Analytical Tendex and Vortex Fields for Perturbative Black Hole Initial Data 
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Tendex and vortex fields, defined by the eigenvectors and eigenvalues of the electric and magnetic 
parts of the Weyl curvature tensor, form the basis of a recently developed approach to visualizing 
spacetime curvature. In particular, this method has been proposed as a tool for interpreting results 
from numerical binary black hole simulations, providing a deeper insight into the physical processes 
governing the merger of black holes and the emission of gravitational radiation. Here we apply this 
approach to approximate but analytical initial data for both single boosted and binary black holes. 
These perturbative data become exact in the limit of small boost or large binary separation. We 
hope that these calculations will provide additional insight into the properties of tendex and vortex 
fields, and will form a useful test for future numerical calculations. 
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I. INTRODUCTION 

The first dynamical simulations of the inspiral and 
merger of binary black holes [IH3] marked a significant 
breakthrough in the field of numerical relativity and 
an important step towards understanding the two-body 
problem in general relativity. They also produced the 
first reliable predictions of the gravitational wave signals 
emitted in these events, which are needed for the analysis 
of data from gravitational wave interferometers. Soon af- 
ter these initial calculations, which adopted equal-mass 
and non-spinning binaries, simulations of binaries with 
unequal masses or non-zero spin revealed important phe- 
nomena, including orbital hang-up (see, e.g., |3]), spin- 
flip (e.g. [5]) and black hole recoil (e.g. [6HI0]: see also 
[TT] for a review). 

In an effort to gain a better physical understanding 
of these phenomena, several researchers have developed 
tools for the visualization and interpretation of gravita- 
tional fields and spacetime dynamics. Here we focus on 
tendex and vortex fields (see [HHH]), which are defined 
in terms of the eigenvectors and eigenvalues of the elec- 
tric and magnetic parts of the Weyl curvature tensor (see 
[TT51 - H7] for an alternative "cross-correlation" approach). 
As we will explain in more detail below, the tendex fields 
describe tidal stretching or compression, while the vortex 
fields describe precession. 

Tendex and vortex fields and their properties have al- 
ready been explored for a number of different types of 
spacetimes. In j!3) . the authors focus on applications to 
weak-field systems and gravitational wave generation by 
such systems, including the case of a Newtonian slow- 
motion binary. Several numerical simulations of binary 
black hole systems are discussed in [T3], where the new 
tools are applied to analyze an "extreme-kick" merger. In 
[14] . the authors use tendex and vortex fields in conjunc- 
tion with results from topology to investigate asymptotic 
properties of gravitational radiation. 

In this paper we add to this list an analysis of ten- 
dex and vortex fields for approximate but analytical ini- 
tial data describing boosted and binary black holes (see 



[15]). The initial data are derived as perturbations of 
Schwarzschild black holes, and become exact in the limit 
of small boost or large binary separation. Given that the 
initial data are analytical, we can also find expressions 
for the tendex and vortex fields, providing an analytical 
example of these fields for a strong-field binary. 

Our paper is organized as follows. In Sec. [n] we briefly 
review the notion of tendex and vortex fields, following 
|12) . and apply these tools to a Schwarzschild black hole. 
We review the perturbative initial data for single boosted 
and binary black holes in Sec. Ill In Sees. IV and [V] we 
apply the new visualization tools to these perturbative 
initial data, relegating some details of the calculations to 
several appendices. We conclude with a brief discussion 



and summary in Sec. VI Throughout this paper we use 
geometrical units with G = c = 1. 



II. TENDEX AND VORTEX FIELDS 

Owen et al. [T2] introduced several tools for visualiz- 
ing spacetime curvature, based on the eigenvectors and 
eigenvalues of the electric and magnetic parts of the Weyl 
curvature tensor. In Sec. II A we follow [T2TTH] and de- 
fine tendex and vortex fields as well as related quantities 
needed for the rest of the paper. In Sec. |IIB| we ap- 
ply these definitions to a Schwarzschild black hole, which 
serves both as a pedagogical example and as the back- 
ground solution in later parts of this paper. 



A. Definitions 

Assuming maximal slicing and vacuum, as we will 
throughout this paper, we can write the electric and mag- 
netic parts of the Weyl curvature tensor as 



Si 



and 



Btj = €j D k K u , 



(1) 



(2) 



2 



respectively. Here lowercase Latin indices run over the 
three spatial coordinates, Rij is the three-dimensional 
Ricci tensor, Kij is the extrinsic curvature, denotes 
the covariant derivative compatible with the spatial met- 
ric 7y, eijk is the spatial Levi-Civita tensor, 



£ijk = l 1,2 [ijk] 



(3) 



where [ijk] is the alternating symbol with [123] = +1, 
and 7 is the determinant of the metric. Maximal slicing 
implies that the trace of the extrinsic curvature vanishes, 



K = ?3Ki. 



0. 



Evaluating the electric part of the Weyl tensor on the 
horizon of a black hole and contracting it twice with the 
horizon's inward unit normal vector N % yields the horizon 
tendicity 



(4) 



The horizon tendicity measures the strength of the tidal 
acceleration at the horizon 12J. The analogous quan- 
tity for the magnetic part of the Weyl tensor yields the 
horizon vorticity 



B 



NN 



BijIPN 3 , 



(5) 



which measures the strength of the frame-drag angular 
acceleration at the horizon [12] , Horizons can be col- 
ored according to the horizon tendicity or vorticity to 
help illustrate their gravitational properties. We provide 
an example in Fig. [I] below, but also refer to the many 
examples in Refs. [12H14| . 

Both £ij and Bij are symmetric, and therefore can 
be characterized by their three orthonormal eigenvectors 
and the associated eigenvalues. It is convenient to work 
in an orthonormal basis where the spatial metric is 5%j. In 
such a frame, which we will denote with hats, we do not 
need to distinguish between contravariant and covariant 
indices. The eigenvectors and eigenvalues satisfy familiar 
equations of the form 



FjV 3 = Xv l 



(6) 



which is equation (3) of [14j except that we raised the free 
index using the spatial metric. We then refer to the eigen- 
vectors of £ij as the tendex fields and the corresponding 
eigenvalues as tendicities, while the corresponding quan- 
tities for Bij are called the vortex fields and vorticities. 
Finally, we refer to the integral curves of the tendex and 
vortex fields as tendex lines and vortex lines. Both £jj 
and Bij are also traceless, so that, in an orthonormal ba- 
sis, their eigenvalues have to add to zero. In the following 
we will derive Bij and B^ in a coordinate basis and then 
switch to an orthonormal basis to solve the eigenvalue 
problem. 

An observer oriented along a tendex vector will tend 
to be tidally stretched for negative tendicity and tidally 
compressed for positive tendicity [TH H3] ■ Neighboring 
gyroscopes oriented along a vortex vector will exhibit 



counterclockwise differential precession for negative vor- 
ticity and clockwise differential precession for positive 
vorticity [TH [T3] . In the following Sections we will repre- 
sent tendex and vortex fields in two-dimensional plots 
with the help of "iron filings" that are familiar from 
representations of magnetic field lines. The iron filings 
show the direction of the eigenvector; we simultaneously 
shade the plot backgrounds according to the correspond- 
ing eigenvalues (see [T^ - TTi] for many examples of alter- 
native representations of these fields). 



B. A Schwarzschild Black Hole 

In this Section we compute tendicities and vorticities 
for a Schwarzschild black hole in isotropic spatial coordi- 
nates, both as a pedagogical example and so that we can 
use the results as the background solution in our later 
perturbative treatment. 

Consider a Schwarzschild black hole with bare mass 
M. . For a Schwarzschild black hole the bare mass is equal 
to the black hole's ADM energy Madm or the irreducible 
mass M\„. In this Section we may therefore replace A4 
with either Madm or M- m . We write the spatial metric 
as 



Jij = i> 7ii> 



(7) 



where ip is the conformal factor and jij is the conformally 
related metric. In isotropic coordinates, on a slice of 
constant Schwarzschild time t, the conformal factor for a 
Schwarzschild spacetime then takes the form 



i/>=l + 



M 

2r ! 



(8) 



and the conformally related metric is flat; in spherical 
polar coordinates we have 



jij = diag(l,r 2 ,r 2 si" 2 ■ 



(9) 



Also, on a slice of constant Schwarzschild time t the ex- 
trinsic curvature vanishes identically, 



= 0. 



(10) 



From ([2]) we see that the magnetic part of the Weyl tensor 
is also zero, 



Bij = 0. 



(11) 



As one might expect for a non-rotating, static black hole, 
all vorticities therefore vanish identically. 

The tendicities of a Schwarzschild black hole, on the 
other hand, are non-zero. In order to evaluate them, we 
first note that for a vanishing extrinsic curvature, the 
electric part the Weyl tensor ([!]) reduces to the Ricci 
tensor 



P. 

^"1 



R, 



(12) 



3 



We can evaluate the Ricci tensor using 

Rij = R ij -2(DiD j lmlj + j ij ^ m D l D m lmlj) 
+ 4((A ln^^-ln^)- 

7 y 7' m (AlnV)(-D m ln^)) (13) 

(see, e.g., [H]), where R\j and Di are the Ricci tensor 
and covariant derivative associated with the conformally 



related metric 



Hi 



Since the latter is flat, we have Rij = 



0, and we find that the only non-zero components of £ij 
are 



£ rr = 



2M 



and 



Be 



sin 2 6 



M 

ip 2 r 



(14) 



(15) 



To compute the horizon tendicity, we first note that 
the inward normal on a sphere of constant radius r is 



V ' =| -^' ' 



Wc then have 
£nn 



■ £ i:j N l N 3 



£ rr N r N r 



2M 



(16) 



(17) 



which, when evaluated at the horizon r — M/2, yields 

1 



£ 



NN 



4M 2 



(18) 



where is the two-dimensional Ricci scalar for the 

horizon and the last equality is a consequence of a more 
general result for quiescent black holes given in Ref . [15] . 

To solve the eigenvalue problem in equation ([6]), we 
convert £u to a spherical polar orthonormal basis. The 
same transformation to an orthonormal basis will be used 
in later sections - only the conformal factor will be dif- 
ferent when we consider perturbed data. Noting that the 
spatial metric is defined by 7y = • a,, we now define 
orthonormal basis vectors 



1 



1 



1 



-eg, 



(19) 



ip 2 ' 6 tp 2 r ' * ip 2 r sin I 
with a corresponding dual basis of 1-forms 

Q f = ifj 2 dr, Q § = xjj 2 rd6, = V> 2 r(sin6»)d0. (20) 
The orthonormal components £%j can be identified using 

£ijdx dx = £ij£j l Cj3 . (21) 

For example, 

1 



£, 



^ V 4 rsin 



(22) 



The other components of are similar, and Bij trans- 
forms just like Sij. Applying these transformations to 
( 14 ) and (fTBl we find that the non-zero components of 



£ij in a spherical polar orthonormal basis are 



2M 



(23) 



This is consistent with the results in [T3] and [T5], but 
expressed in a different coordinate system. Equation ( 23 ) 



shows that £%j is diagonal in this basis, so its orthonormal 
eigenvectors can be chosen to be 



,(o) * 



= (e r f = (1,0,0), 



and 



4f = (e e -)*= (0,1,0), 



= (*£)' = (o,o,i). 



(24) 
(25) 

(26) 



Here the superscript has been added for consistency 
with later Sections where these results will serve as back- 
ground solutions to a perturbative treatment, and the 
labels El, E2 and E3 denote these as the three eigen- 
vectors of £%j. The corresponding eigenvalues are given 
by 



A 



(0) _ £_(0) _ 



El 



2M 



and 



\(0) _ X (0) _ o(0) _ 
A E2 — A E3 — °ao ~ 



Since A^j is negative, observers are stretched in the radial 

direction. The eigenvalues A^ an d A^?] are positive and 
degenerate, so observers are compressed equally in all 
tangential directions (see also Fig. 2 in cither of Refs. 
[T2l [T3l and the associated discussion). 



(o)' 



M 



128Mr 3 
(M + 2r) e 



6-iMr 3 



tfof 3 {M + 2rf 



(27) 



(28) 



III. PERTURBATIVE BLACK HOLE INITIAL 
DATA 

Numerical relativity simulations using a "3-1-1" decom- 
position require initial data, namely a spatial metric 7™ 
and extrinsic curvature Kij satisfying the Hamiltonian 
constraint and the momentum constraint. Using the con- 
formal transformation Q , and transforming the extrinsic 
curvature according to 



Kij — i)j Aij, 



the Hamiltonian constraint takes the form 

1 



D 2 Tp = --Tp- 7 A i jA i i, 



(29) 



(30) 



4 



while the momentum constraint reduces to 



Here 



0. 



(31) 



Here Di denotes the covariant derivative associated with 
the conformally related metric 7y, and D 2 = ■j lJ DiDj is 
the associated Laplace operator. In the above expressions 
we have again assumed maximal slicing and vacuum; we 
have also assumed conformal flatness, meaning that "fij 
is a flat metric and D 2 a flat-space Laplace operator. 

The momentum constraint has become linear under 
these assumptions, and is solved analytically by so-called 
"Bowen-York" solutions [2"Ul422| . For a black hole at co- 
ordinate location C l with linear momentum P l these so- 
lutions are given by 



2r% 



pV c 



(f- 



i(,)P k n k c 



(32) 

where tq — \x l — C*| is the coordinate distance from the 
center of the black hole and n l c — (x l — C l )/rc is the 
unit vector (normalized with respect to the conformally- 
related background metric) which points from the center 
to coordinate location x % . 

The solutions ( |32|) ca n then be inserted into the Hamil- 
tonian constraint (30), which, in general, still has to 

for different ap- 



be solved numerically (see, e.g., 
proaches and results, as well as [TT] for a review.) Here 
we will review a perturbative but analytical approach 
(see [T51 HS1 ISO])- We note that for vanishing boost, 
the Hamiltonian constraint ( 30 1 is solved exactly by the 
Schwarzschild conformal factor For non-zero boost, 
we can then consider the leading-order perturbations of 
the Schwarzschild conformal factor. In the following two 
Sections we will consider single boosted black holes and 
binary black holes separately. 



A. Single Boosted Black Holes 

We first note that A*- 7 ' is linear in the magnitude P of 
the momentum. From the Hamiltonian constraint (30) 



we then see that all perturbations of ip can only contain 
even powers of P. Defining 

P 

€P = M> 

we can therefore write the solution to the Hamiltonian 



(33) 



constraint (30) as 



(34) 



where V'(o) is the Schwarzschild conformal factor (|8j) for 
a black hole at coordinate location C\ 

V>(0) = 1 + ^, (35) 

and where the function u can be written as 

M ( 

u {r c )Po{cos9) +u 2 {r c )P2{cos6) 



8(M + 2r c 



and 



P (cos6>) = 1 



P 2 (cos6») = 2( cos2 1 



(37) 



(38) 



are Legendre polynomials with 9 measured from the 
boost direction, and the radial functions Mo(^c) an d 
u 2( r c) appearing in equation (36) are given by 



u (r c ) =M 4 + lOM 3 r c + ^QM 2 r 2 c + 80Mr^ + 80r£, 

(39) 

and 



U2(r c ) 



., , 42X 5 r c + 378X 4 r^ + 1316A^ 3 r^ 
+ 1\hZM 2 r% + 1536Xr^, + 240r£ 
+ ^M(M + 2r c fln( j4 ^ 2 -) \ (40) 



(see [291 [30]). 



B. Binary Black Hole Systems 



The initial data discussed in Sec. MI Al can be used 
to construct similar initial data for a binary black hole 
system with one hole at coordinate location C\ with bare 
mass M i and momentum PJ , and another at coordinate 
location C\ with bare mass M.i and momentum P| (see 



The momentum constraint (31) is linear in A lJ , so it 



can be solved by adding two solutions of the form shown 
in equation ( 32 ) : 



71 — ^CiPi 



SA P . 



(41) 



The Hamiltonian constraint (30) is not linear in ip, but 



it is still possible to construct a perturbative solution 
similar to that for a single boosted black hole by adding 
the separate perturbations from each of the holes. The 
derivation in |18j assumes that the coordinate distance s 
between the holes obeys 



s > Mi > Pi, 
s > M 2 > P2, 



(42) 



with Pi/Mi ~ P2/M2 ~ ep- The conformal factor sat- 
isfying equation ( 30 ) to second order in the boosts is 



^ = l + ^ + ^+e 2 P ui + e 2 P u 2 +0(e P ), (43) 
where Ui and u 2 are appropriate versions of equation 



(36) (36), with rc replaced by either rc x or rc 2 



5 



In this paper we will assume an equal-mass binary with 
Mi = M.2 = M.. We also assume a quasicircular orbit, 
which implies that the black holes have equal magnitude 
momenta P\ = P2 = P, so ep 1 = ep 2 = ep. Applying 
the virial relationship from equation (57) of [TS] results 
in the Kepler law 



M 

2s 



(44) 



Note that in using equation (57) of [18] we have replaced 
the irreducible mass |31| with the bare mass, but this 
difference only changes the calculation at order O(ep), 
so it can be ignored here. 

For both single boosted and binary black holes, A^ is 
linear in P, so that perturbations of ip contain only even 
powers of P. From (29) we then see that Ky can only 
have odd powers of P. Therefore, perturbations of the 
Schwarzschild expressions for £ « also contain only even 
powers of P, while perturbations of Bij contain only odd 
powers. 



IV. TENDEX AND VORTEX FIELDS FOR 
BOOSTED BLACK HOLES 

We now consider a single black hole at the origin 
boosted in the +z direction, described by the pertur- 
bative initial data from Sec. |III A| We then compute 
the tendex and vortex fields to leading-order, adopting 
the results of Sec. 



solution. 



II B as the unperturbed background 



and 



'(2) 



4« 



,(0) 



ip^r 2 sin 1 



^(o) r2 s ^ n2 ' 



(49) 



where the function u is given by equation ( 36 1 



B. The Magnetic Part of the Weyl Tensor 

Only odd powers of ep appear in the magnetic part of 
the Weyl tensor, so, to leading order, it can be written 
in the form 



Bm = e P B$ + 0(ep). 



(50) 



In Appendix |B 2| we show that the only non-vanishing 
components of the leading-order term By , expressed in 
a coordinate basis, are the rcj) and <f>r components 



„(i) _ „(i) _ 967Wr 3 sin 2 6 
r * ^ ~ (M + 2r) 5 ' 



(51) 



We then use the transformation p2| to find the orthonor- 
mal components 



(1) _ r(1) _ 



B\ L = B 



1 



' f -0(o) rsm 



:B 



(1) 



(52) 



C. Horizon Tendicity and Vorticity 



A. The Electric Part of the Weyl Tensor 

As we have discussed above, only even powers of ep 
can appear in perturbations of the electric part of the 
Weyl tensor, so that we may write the result as 



c _ 0(0) 



e P £^+0(e%). 



(45) 



For a black hole at the origin we can replace rc with 
and find that in a spherical polar orthonormal basis, £ ~~ 



(0) 



has nonzero components given by equation ( 23 ) with the 



conformal factor ip replaced by V'(o)- The perturbations 

(2) 

£y , expressed in a coordinate basis, are computed in 



Appendix |B 1| (see equations ( B8||B11 )). We then obtain 
the perturbations of the orthonormal components from 



£ 



(2) 



Vfo) 



£(2) 



-ill 



p(o) 



(46) 



(0) 



In Appendix B 3 we adopt the results for £ij in a co- 



ordinate basis to find the horizon tendicity. This cal- 
culation also relies on results for the perturbed horizon 
location that were found in [18j. We compute the hori- 
zon tendicity directly, following the procedure outlined 
in Sec. |II A| but in an alternative approach it is also pos- 
sible to employ the Newman-Penrose formalism (see 
|13|). The result of this calculation is 



1 



1 



■*NN 



16M 2 



(53) 



(-1871 + 2688 In 2) 
' 64CLM 2 



P 2 {cos9)+0(e P ). 



The first two terms simplify if we express results in terms 
of the black hole's irreducible mass. From equation (26) 
of [18] we have 1 



M = M irr 1 



P 2 



8Md 



0(e%) 



(54) 



c(2) 



c(2) _ 



?(2) 
J r9 ' 



7(2) 



(0)' 



4u c(0) 

5 r 2°ee > 
(o) r 



(47) 



(48) 



1 We note that this result holds only for wormhole data. For trum- 
pet data, for example, the bare mass M appears to be equal 



to the irreducible mass 

EH El. 



for single boosted black holes; see 



6 



M 



-as 




is 



nS 



FIG. 1: The deviation from the average horizon tendicity for 
a single black hole boosted with P = QAM in the positive z 
direction (i.e. pointing up). The horizon tendicity is negative 
everywhere, but it is more negative (darker shading) near the 
poles at 6 = and 8 = n and less negative (lighter shading) 
near the equator at 9 = n/2. The true distorted shape of the 
horizon is shown, but for the small boosts relevant for this 
paper the distortion from the unperturbed spherical shape is 
not readily apparent. 



Inserting this relation into the expression above we see 
that the (proper area weighted) average value of the hori- 
zon tendicity is 



4M? rr 



o(4). 



(55) 



The simplicity of this result is not surprising on physi- 
cal grounds and we might have anticipated it from the 
Gauss-Bonnet theorem as discussed in this context in 
[12]. As shown in Appendix |B 3[ the average horizon 
tendicity for a horizon with spherical topology will take 
the above form whenever a particular combination of spin 
coefficients vanishes when integrated over the horizon - 
which is true here to the relevant order. 

To leading order in ep, the deviation from this average 
value is proportional to P2(cos8). In Fig. [I] we plot this 
deviation of the horizon tendicity from its average value 
for a black hole boosted in the +z direction with P = 
0.1M. 

In contrast to the horizon tendicity, the horizon vor- 
ticity vanishes, at least to the order of our analysis. Fol- 
lowing the same steps that were used to derive £nn m 



b™ = b« 7 £ ) + 0(4) = 0(4), 



(56) 



because from equation ( B13 1 the only nonzero compo- 



nents of Bij to this order are B r ^ and B^ r . 



D. Tendex Fields 

In order to find the tendex fields we have to find the 
eigenvectors and eigenvalues of the electric part of the 
Weyl tensor. We have already solved this problem for the 



unboosted background solution in Sec. II B We can now 
find the leading-order corrections to the Schwarzschild 
results by solving the eigenvalue problem perturbatively. 
While this "stationary perturbation theory" technique is 
well known, especially in the context of quantum me- 
chanics, we summarize the most important results, and 
specialize to the three-dimensional matrices encountered 
in this context, in Appendix [XJ Given that two of the un- 
perturbed eigenvalues are degenerate (see equation (28)), 
both degenerate and nondegenerate perturbation theory 
are required. 

Applying these techniques to the perturbations of the 
electric part of the Weyl tensor as derived in Sec. |IV A| 
above, we find that the nondegenerate eigenvalue (27) 
generalizes to 



(57) 



Here the background term \g\ is given by equation ( 27 ) 
and 



X (2) _ 



4M 



2LW 2 (3(cos2(9) 



5(M +2r) 12 

x (.M 2 + 2Mr + 24r 2 ) | M + 2r 

2TW 7 + 231M 6 r + 1540.M V + 69307W 4 r 3 

- 187207W 3 r 4 + 27568.MV + 188167Wr 6 

- 3360r 7 ) - 2r(21A4 8 + 231A4 7 r + 1550.MV 



693CLM 5 r 3 + 18121XM V + 24368A4V 
11616A4 2 r 6 - 4320A4r 7 + 3200r 8N ) ). (58) 



(59) 



The corresponding eigenvector is 

V E1— V E1 + e P v El + u ( e P)i 

where the background term v^l % is given by equation 



equations (BI5) - (B19), we see that the horizon vorticity (24) and where the only non- vanishing component of 



7 



(2) t . 
«B1 1S 



,(2)0 



<E1 



80r 3 (X + 2r) 6 



3840r s sin 6 cos ( 



sin 20 ) LM + 2r ) 2r{ 21M 6 + 357M 



2170M 4 r 2 + 6342X 3 r 3 + 93887W 2 r 4 
62l6Mr 5 + 720r 6 ) +21.V|(A1 + Nr 

5 / X 



(.M + 2r) ln(- 



(60) 



The degenerate unperturbed eigenvalues become 

XE2 = \ { E ] 2 + € 2 P \f 2 +0(e%) (61) 



and 



(62) 



where the identical background terms A^ an( i A^g are 



given by equation ( 28 ) . The two perturbative corrections 
are 



A 



(2) 
E2 



AM 



and 



Al = 5(A , + 2r)12 QMC+2TD), (64) 

where the coefficients A, B, C and D are given by 

- M \ 

(?>( cos 26») (M 2 + 3Mr + 14r 2 ) 

- (M-r)(M+6r)^j 

- 2r( cos 2fl) (2LM 7 + 252M 6 r + 1519A4 V 

- 5698.MV + 132167W 3 r 4 + 17536A4 V 

- 11184.Mr 6 + 2400r 7 Y (65) 



5(M + 2r 



AM 



12 (3MA - 2rB) 



(63) 



A = 7M[M+2r) In 



S = 2LM 8 + 29AM 7 r + U72M 6 r 2 + 3234.M 5 r 3 
+ 2508A4 4 r 4 - 928.MV - 4807W 2 r 6 
+ 4320A4r 7 - 1600r 8 , (66) 

c = 7M { M+2r T h <M^) 

x (2M 2 - 3r ( cos 29) (M - lOr) + 7 Mr + 18r 2 ) 

+ 2r 2 ( cos 29) ( - 21M 6 + 2\M b r + 1232.MV 

+ 550AM 3 r 3 + 10032A4 2 r 4 + 76327Wr 5 

+ 960r 6 ) , (67) 



and 



D = A2M 8 + 52HM 7 r + 3022A4V + 10164AtV 
+ 206287WV + 23440A4 3 r 5 + 11136A4 2 r 6 



8 



+ 1600r 

The corresponding eigenvectors are 



and 



V E2~ V E2 + e P v E2 + u ( e P> 



v E 3=v ( Er+v(4) 



where the only nonvanishing component of v£l 1 is 



(2)f 
' } E2 



(2) e 

' J E1 ■ 



(68) 
(69) 
(70) 

(71) 



The eigenvectors v l E1 , v E2 and v E3 are orthonormal to 
second order in ep as expected. For small boosts, these 
eigenvectors are not visibly different from the purely ra- 
dial or purely tangential unperturbed eigenvectors. Alge- 
braically, we see that the perturbation mixes the and 
eg eigenvectors of the unperturbed state but leaves the 
eigenvector unchanged. This is expected because the 
boost breaks spherical symmetry but preserves axisym- 
metry. 



E. Vortex Fields 

As we have seen above, only odd powers of P enter 
into an expansion of the magnetic part of the Weyl ten- 
sor. For Schwarzschild we have Bij = 0, so even for small 
boosts we cannot treat this problem perturbatively. To 
leading order, Bi-j takes the simple form given by equa- 
tions (50)-|52|), so that it is straightforward to compute 



the eigenvalues and eigenvectors directly by diagonaliza- 
tion. 

We find that the first eigenvalue vanishes to our order 
of analysis, 



A 



Bl 



0(e 3 



pi- 



rn 



while the two other eigenvalues take the values 

1536Xr 6 sin6» m . 3 , 
52 = - X B3 = -e P {M + 2r)9 + 0(e 3 P ), (73) 

The corresponding orthonormal eigenvectors are 

*4i = (e e y, (74) 



and 



J B2 



J B3 



1 
1 



(<*)' + (e 3 )' +0(e P ), 



(<*)'" -(e J*" )+0(e 2 p) 



(75) 



(76) 



Fig. [2] shows Xb2 in the z — plane (orthogonal to the 
boost of the black hole), along with the projection of v l B2 
into that plane. In Fig. [3] we show the same quantities, 
but in the y — plane. 
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FIG. 2: The vortex eigenvalue Xb2 in the 2 = plane for 
a black hole boosted in the positive z direction, along with 
the projection of eigenvector vg 2 into that plane. Randomly 
placed line segments of fixed length mimic the appearance 
of "iron filings" and trace the orientation of the eigenvector, 
while the colored background shading shows the behavior of 
the eigenvalue. The eigenvalue is most negative where the 
shading is darkest and approaches zero as the shading fades 
away. In this figure the black hole is boosted in the direction 
pointed out of the paper, at the reader. Note that this figure 
and later figures bear some resemblance to the ridge patterns 
in Fig. 3 of [HI. 



TENDEX AND VORTEX FIELDS FOR 
BINARY BLACK HOLES 



We now turn to perturbative initial data describing 
an equal-mass binary black hole system in quasicircular 
orbit. As demonstrated in [15], these solutions can be 
constructed as perturbations of two Schwarzschild black 
holes, by adding the perturbations created by the black 
holes' boosts to those created by the presence of the bi- 
nary companion. Here we follow a similar approach to 
construct the tendex and vortex fields for such a binary 
system. While the calculations can be carried out an- 
alytically at least in principle, some of the expressions 
become very unwieldy and do not provide much insight. 
In some cases we therefore restrict the analysis to certain 
regions or symmetry planes that allow a direct compari- 
son with the results for single boosted black holes in Sec. 

El 



FIG. 3: The vortex eigenvalue As2 in the y = plane along 
with the projection of eigenvector v l B2 into that plane. In this 
figure the boost of the black hole is pointing upwards. 



A. The Electric and Magnetic Parts of the Weyl 
Tensor 

The perturbations of the electric part of the Weyl ten- 
sor can be constructed as for a single boosted black hole, 
adding contributions from the two black holes' boosts and 
then taking into account the perturbation created by the 
binary companion. The resulting expressions, however, 
are very messy, and we therefore restricted this analysis 
to regions either close to one of the black holes, or far 
from both. We will return to this analysis later. 

To compute the magnetic part of the Weyl tensor we 
can also add the perturbations created by the two holes. 
As shown in Appendix |C 1[ it is convenient to first con- 
vert the results for a single boosted black hole to a Carte- 
sian orthonormal basis with the origin placed at the cen- 
ter of mass. In this coordinate system, the two black 
holes, separated by a coordinate distance s, are placed at 
coordinate locations (±s/2,0, 0) with boosts (0, ±P, 0). 
Keeping symmetry in mind, the nonzero components of 
By, to the desired order in P, are 



E- - 



2z 

y 

A 



r 4 Cl (s-2x) 



e P 1536Mz x 



(77) 



r%(s + 2x) 



{M + 2r Cl f (M+2r C2 y 



+ 0(e P ), 



= epl536M.yz 



(M + 2r C2 ) 9 (M + 2r Cl f 



0(e 



(78) 
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and 



r 4 c (s-2x 2 -4z 2 

tp384M Cl y. . ' , Q 

I {M + 2r Gl f 



Az 2 



(M + 2rc 2 ) 9 



+ 0(e%). (79) 



B. Horizon Tendicity and Vorticity 

To calculate the horizon vorticity it is convenient to 
change coordinates and consider Bij in the vicinity of one 
of the holes (say hole 1 at coordinate location (s/2, 0, 0)). 
In a spherical polar orthonormal basis centered on this 
hole, with the boost in the (new) z direction, the leading- 
order components of Bij are just those specified for the 
single boosted black hole in equations ( 50 1-( 52 ) - this is 



because Kepler's law (44 1 implies that l/si scales with e|>, 
so that the corrections created by the companion black 
hole are of higher order. As shown in [18], corrections to 
the horizon location due to the companion black hole are 
also of higher order, so the horizon is still axisymmetric. 
As a consequence, the horizon vorticity still vanishes to 
our order of analysis. 

To find the horizon tendicity for hole 1, say, we first 
need to find expressions for the components of £ij near 
hole 1 - this is done in Appendix C 2 We expand £g as 



we did for the single boosted black hole in equation ( 45 1 . 



As before, the background is still given by equation 



(2) 

(23) with ip replaced by t/j^y The perturbations £^' 
expressed in a coordinate basis are given by equations 
( C8|C11 ) in Appendix | C 2| - they are very similar to those 
found for a single boosted black hole, but differ in some 
terms because of the presence of the companion black 
hole. 

From £ij we compute the horizon tendicity following 
the calculation in Appendix [B] for single black holes, ex- 
pect that we replace the conformal factor with that for a 
binary. The result is 



£ 



1 



9 



NN 



4M : 
e 2 ( 



^16M 2 
-1871 + 2688 In 2) 



(80) 



640A4 2 



P 2 (cos0) + 0(4). 



As for the single boosted black hole, we can simplify the 
first two terms by expressing them in terms of the irre- 
ducible mass. The relationship between the bare mass 
and the irreducible mass in a binary is given by equation 
(47) in [IE], 



M = M irr 



M 2 

^+0(4). (81) 



8M irr 

We therefore see that the average horizon tendicity is 
again 

1 



cave 
C NN 



AM 2 



+ 0(4). 



(82) 



^v+, •sj 



/! 



/ 



-__ ;+ 



•> — , ■ 




o 



FIG. 4: The vortex eigenvalue \bi in the symmetry plane 
x — between the two binary companions, along with the 
projection of eigenvector v' B2 into that plane for a binary at 
separation s — 50 M. The analogous figure for Xb3 and v B3 
would look like this figure flipped vertically about z — 0, 
but with the degree of shading indicating the magnitude of a 
positive eigenvalue. 



To leading order, the deviation from this average value 
is identical to that for a single boosted black hole, so a 
plot of the deviation would look like Fig. [T] The spin 
coefficient argument made in Appendix |B 3| for the single 
boosted black hole applies here as well - equations (B24 



B30 ) are unchanged even though expressions for their 
constituents change slightly. 



C. Vortex Fields 

While the eigenvalues and eigenvectors of the magnetic 
part of the Weyl tensor can be computed everywhere, 
they are in general quite complicated and do not offer 
much insight. We instead restrict our analysis to certain 
planes, some of which allow for a direct comparison with 
the results for a single boosted black hole. 



1. The symmetry plane between the black holes 

The x = symmetry plane midway between the black 
holes provides the simplest case. Defining tq = rc t — 
rc 2 = \J (s/2) 2 + y 2 + z 2 we can write the eigenvalues as 

m2Mr% sz 
B1 ~ P [M + 2r c ) 9 ' ( P >> ' ( ' 
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1536A44sMVg+g) , 2 
A b 2 -^p ( ^ + 2rc)9 + 



(84) 



and 



AB3 = €P (M+2r c r (1 + ° {ep)) ■ 



The corresponding eigenvectors are 

v% 1 = (1,0,0) + o(4,), 



J B2 



(85) 
(86) 

(87) 



and 



v i B3 = A(o,-(z+Vy Y T^),y^ + 0(e 2 P ), 
where the normalization factor A is given by 



.4 



V2^y 2 + z(z + vV + z^) 



(89) 



In Fig. [4] we show the eigenvalue Amhithea; = sym- 
metry plane along with the projection of eigenvector v B2 
into that plane. In this figure, and in all the following 
figures, we assume a binary separation s — 50.M. 



-20 



o 
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FIG. 5: The vortex eigenvalue As2 in the symmetry plane 
y = orthogonal to the black holes' boosts, along with the 
projection of eigenvector Vg 2 into that plane. A close-up view 
of the region near hole 1 is shown in Fig. [6] 



where the normalization factor A is given by 



2. The symmetry plane orthogonal to the black holes ' boosts 

Expressions for the eigenvalues and eigenvectors in the 
y = and z — planes are somewhat more complicated 
but can be written compactly in terms of components of 
Bij. For the y = symmetry plane, which is orthogonal 
to the black holes' boosts, the eigenvalues are 



Am = 0(4), 



\ B 2 = -V(B*x) 2 + (Bxi) 2 + 0(4), 



and 



a B3 = sfW&f 2 + (£ii) 2 + 0(4)- 

The corresponding eigenvectors are 

v% 1 = (0,1,0) + 0(4), 

v% 2 = A(B ££ - V / SF+W,0,Bi i ) 



+ 0(4), 



and 



y B3 



(90) 
(91) 

(92) 

(93) 

(94) 



+ 0(4), (95) 



V2j(B ££ )i + {B & ,y + B ££ ^{B~ £ 



(96) 

In Fig. [5] we show the eigenvalue A#2 hr the y — plane, 
together with the projection of eigenvector into that 
plane. A close-up view of the region near hole 1 is shown 
in Fig. [6] 



3. The orbital plane 



Finally, in the orbital plane z = the eigenvalues are 



Am = 0(4) 



\B2 



7(£x*) 2 + (%) 2 + 0(4), 



and 



A B3 = v /(6 fii ) 2 + (Kp) 2 + 0(e? 3 ) 
The corresponding eigenvectors are 

v* B1 = AV2(~By i ,B £i ,0)+O(e 



v B2 = A J (B«) 2 + (B fl< ) 



0(e 



(97) 
(98) 

(99) 

(100) 

(101) 
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FIG. 6: Same as Fig. [5] but in the vicinity of hole 1. Compare 
to Fig. [2]- note that hole 1 is boosted into the page here while 
the hole in Fig. [2] is boosted out of the page. 



FIG. 7: The vortex eigenvalue Xb2 in the orbital plane z = 
along with the projection of eigenvector Vg 2 into that plane. 
A close-up view of the region near hole 1 is shown in Fig. [8] 



and 



v 



{Byz? 



where the normalization factor A is given by 

1 

2, 



.4 



^xz J 



(Byz) 2 



(102) 



(103) 



Fig. [7] shows eigenvalue A 32 in the plane z = along 
with the projection of eigenvector v l B2 into that plane. A 
close-up view of the region near hole 1 is shown in Fig. 

m 



D. Tendex Fields 

While the electric part of the Weyl tensor can, in prin- 
ciple, be computed everywhere, it is in general quite com- 
plicated, meaning that the expressions for its eigenvalues 
and eigenvectors are even more unwieldy. We therefore 
restrict our analysis to two regions, namely the vicinity of 
one of the black holes, and the asymptotically far region. 
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Same as Fig. [7] but in the vicinity of hole 1. Compare 



E 



1. Vicinity of one black hole 

As before, we use the perturbation theory results sum- 
marized in Appendix [A] to find the perturbed eigenval- 
ues and eigenvectors of near hole 1. We change co- 
ordinates as in the horizon vorticity calculation above, 



replacing rc 1 with r for notational convenience. To sec- 
ond order in ep, the perturbed eigenvectors are exactly 
those given by equations ( 59 1 , ( 69 ) and ( 70 ) for the single 



boosted black hole, but the eigenvalues are slightly dif- 
ferent for this problem. The expansions (57 1, (61), and 



(62 1 apply as before, but the perturbative pieces change 
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slightly. The perturbative piece of the first eigenvalue is 
given by 



x (2) 

A E1 



AM 



5(M + 2r) u 



21M 



' (3 (cos 26») 



1 



x ( M 2 + 2Mr + 24r 2 ) ( M + 2r 
M 



x In 



6Mr{ cos 26) x 



kM + 2r 

(2.IM 7 + 231M 6 r + 1540AtV + 6930MV 

+ 18720.M 3 r 4 + 27568A4V + 18816A4?' 6 
+ 3360r 7 ) +2r( - 21M 8 - 23LW 7 Y 



- 1630XV - 693(XM 5 r 3 - 13320A4 V 
+ 1232X 3 r 5 + 45984A4 2 r 6 + 657607Wr 7 

+ 22400r 8 ^. (104) 

The perturbative pieces of the second and third eigen- 
values are given by equations ( 63 ) and ( 64 ) . The auxil- 



iary variables A and C remain the same as in the single 
boosted black hole case, but B and D are now given by 



B 



21M & + 294M 7 r + 14327W 6 r 2 
+ 4908.M 4 r 4 + 11872.M 3 r 5 
+ 35040A4r 7 + 11200r 8 



5„3 



3234A4 r 
28320\A/l 2 r 6 



and 



D 



42M 8 + 5257W' 
+ 18228.MV 4 



r + 3062X 6 7 
10640M 3 r 5 



10164MV 



(105) 



5„3 



17664.M 2 



30720Xr 7 - 11200r 8 



2. Asymptotic region 



(106) 



We can also find the eigenvalues and eigenvectors of 
£ij for the binary far from both holes. We return to 
our original coordinate system for the binary, with r = 
\/ x 1 + y 2 + z 2 measured from the center of mass. For 



r > s, we can expand the conformal factor (43 1 in powers 
of e r = s/r to find 



1 



,5^ 



^ = l + -\2M + 2"—j+0{e%e r ) (107) 

(where P is the magnitude of each individual black hole's 
momentum). We can therefore identify the ADM mass 
of this system as 



Madm =2[M 



5P 2 \ 
8M ) 



•0(4) 



(108) 



(see also equation (44) in [E]), where the first term on 
the right hand side accounts for the bare masses of the 
two black holes, and the second term for the binding 



and kinetic energies. To leading order, therefore, the 
conformal factor of the binary appears like the conformal 
factor of a single black hole with the above ADM mass 
located at the origin. 

To compute the electric part of the Weyl tensor, we 
recognize that the leading order terms in the Ricci ten- 
sor Rij arise from the ADM mass term in the conformal 
factor above, and that the contributions from the extrin- 
sic curvature fall off more rapidly and can be neglected. 
To leading order, therefore, the electric part of the Weyl 
tensor is identical to that of a single Schwarzschild black 
hole, if we express its mass in terms of the ADM mass. 



Borrowing the results from Section II B we have 



2 Ma 



dm 



0(e 



P c r j 



(109) 



Similarly, the tendex fields for the binary in the asymp- 
totic region are identical to those for a single black hole 
derived in Section II B[ just with the bare mass replaced 
by the ADM mass. 



VI. SUMMARY 

The authors of [T^] recently introduced tendex and vor- 
tex fields, defined in terms of the eigenvectors and eigen- 
values of the electric and magnetic parts of the Weyl 
tensor, as an aid to visualize spacetime curvature. In 
particular, the method has promise to help interpret re- 
sults from numerical relativity simulations, and to pro- 
vide insight into the physical processes governing the co- 
alescence and merger of binary black holes as well as the 
emission of gravitational radiation. 

Many numerical simulations of binary black holes start 
with wormholc initial data constructed with the puncture 
method (see, e.g., [27]). Here we present perturba- 
tive but analytical expressions for the tendex and vortex 
fields, based on the perturbative treatment of these initial 
data presented in [18]. In this approach, the boosts of the 
individual black holes and the effect of the binary com- 
panion are treated as perturbations of a Schwarzschild 
black hole, so that the data, and hence our expressions 
for the tendex and vortex fields, become exact in the limit 
of vanishing boost or large binary separation. 

Our results complement other examples of tendex and 
vortex fields presented in [12TU4| . and help us to bet- 
ter understand the properties of these fields. We hope 
that, as analytical expressions for strong-field objects, 
our results may also be useful for comparison with future 
numerical calculations. 
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Appendix A: Finding Eigenvalues and Eigenvectors 
Pert ur bat ively 

Finding the eigenvalues and eigenvectors of a per- 
turbed matrix is a common problem in quantum mechan- 
ics that in its simplest form can be handled by stationary 
perturbation theory. It is less familiar in this context, so 
in this appendix we review well-known results (see, e.g., 
[551 [35]) just to the extent needed for the calculations 
in Appendix [B] and Appendix [C] Note that this is first 
order perturbation theory, but that for consistency with 
the application and notation elsewhere in the paper, we 
work to first order in a small quantity called e 2 . We also 
choose to work in an orthonormal basis. 

Let £.? be a 3 x 3 symmetric matrix with real entries. 



£jj has three real eigenvalues \ E {, A^j an d A^, and 
corresponding orthonormal eigenvectors v E { l , v E l l , and 



(0) x(0) 



,(0) 



Mi 

J E3 



We now consider a perturbation of £~9 of the form 



p.. 



2c(2) 



e z £ 



I] I] 

We then want to find the new eigenvalues 



(Al) 



Aei 


X (0) 
- A E1 ~ 


he 2 Ai- 






Ae2 


_ \(0) 
- A E2 ~ 


He 2 Ai- 


F0(e 4 ), and 




Ae3 


_ \(0) 
- A E3 ' 


He 2 Ai- 




(A2) 



along with corresponding new eigenvectors normalized to 
order e 2 



v 



El 


- „(°) 1 4. P 2„(2) i 

— v E1 +e v E1 - 


^( £ 4 ), 




E2 


- „(°) 1 4. A,(2) i 
— w £2 + e V B2 


hO(e 4 ), and 




E3 


- „(°) J 4. A,(2) S 


hO(£ 4 ). 


(A3) 



Exactly how we do this depends upon the degeneracy of 
the eigenvalues. 



2. Degenerate Eigenvalues 



Suppose that the unperturbed matrix £~9 has a two- 



fold degeneracy, e.g. A^ = A E3- ^ n this case the normal- 
ized vectors v E \ 1 and v E \ 1 might be any two orthonormal 
vectors spanning the corresponding subspace. Choose 
^ 1 and Vill 1 so that e 2 £-a vhl l v^l 3 = 0. This ensures 



(0) 



J E2 



E3 



■'ij U E2 V E3 



that the unperturbed states are indeed the limit of the 
perturbed states as e 2 — > 0. Then, the perturbed eigen- 
values are 



Ab2 — A 



(o) 

E2 



■<?£g>v& t v% 3 + Otf), (A6) 



and 



A 



- \W i c 2p(2) (0)» (0)5 n , <K /a 7 n 

E3 - A B3 + e ^3 u B3 + U(e ), (A7j 



(0) 



and if the perturbation succeeds in breaking the degen- 
eracy at first order in e 2 , the eigenvectors are 



J E2 



,(0) I 
' E2 



g(2) (0)5 (0)fe 



A 



(o) 

£2 



A 



(0) 
El 



(A8) 



and 



u iS3 



„(0) I 
J _B3 



f (2) (0)5 (0)fe 



A 



(o) 

£3 



A 



(0) 
£1 



(A9) 



Appendix B: Single Boosted Black Holes 
1. The Electric Part of the Weyl Tensor 



From the form of the initial data specified in Sec. |III A| 
we can expect the lowest-order perturbation to £y to be 



of order e 2 P : 



_ p(o) 



f- ■ — f y 



4^ + 0(4). 

If we write the Ricci tensor as 



(Bl) 



1. Nondegenerate Eigenvalues 



Ri 



(B2) 



Suppose that X E1 is a nondegenerate eigenvalue. The 

(A4) 



corresponding perturbed eigenvalue is then 

\ - \( 0) ^ c 2c(2)„,(on (0)5 , n / 4\ 
Abi - A E1 + e t l0 v E1 v E1 + U(e j, 

and the corresponding perturbed eigenvector is 



-'El 



,,(°) 5 
El 



A2) (0)3 (0)k 
2 C 3 k V E1 V E2 (0)i 

' ,(0) A (0) V E2 



A 



El 



f (2) (0)5 (0) k 
2 jk El E3 

,(0)_ A (0) V E3 



E2 

(o)i 



X 



El 



+ 0(e 4 ). (A5) 



£3 



we see that 



which we already computed in Sec. |IIB| and 

^=R^-^K^ 



(2) 

The perturbative term £ ^ can then be found by substi- 
tuting equation (34 1 into equations ([7]), ( p9| ), and (13 1 
We note that 



(B3) 



(B4) 



7(o) 



(B5) 



14 



K W - ? /,- 2 z.. 



(B6) 



and wc compute the perturbation of the Ricci tensor from 



u 



■>(n,or • 



ln^(o))D. 



■j(^( } u ) 



+ 4f v id) ij-'j > ' , 

+ Dj (In ■i/j i0) )D i (if>( Q 1 ) u) 

+ D m (lnV( ))A(^ju))). 



(B7) 



of £ i 2) to be 



£(2) 



From ( 36 ) and ( 32 1 we then find the nonzero components 



M 



2Qr i (M+2rf 



3M( - 2r(cos29)(21M 7 



+ 399M 6 r + 30527W V + 12194.M V 
+ 27344MV + 33712.M 2 r 5 + 19776Mr 6 
+ 3360r 7 ) - 7M (3(cos 29) +1)(m + 4r 

M 



(m + 6r) (M + 2rj In ( 



M + 2r 



- 2r ^2LM 8 + 399M 7 r + 3062M e r 2 

+ 12274MV + 275UM 4 r 4 + 33712X 3 r 5 

+ 189767W 2 r 6 + 208OMr 7 + 3200r 8 ) ) , (B8) 



r(2) _ 
'66 ~ 



M 



20r 2 (M + 2r) s 
M 



3M (jM(M + 2 r y 



7Mr 



5 r 2 



and 



,(2) 



+ 14r 2 ) - (M - 2r) (M + 3r)) 

+ 2r(cos 29) (21M 7 + 336.M 6 r + 2275.M 
+ 8330MV + 17528X 3 r 4 + 20608.M V 
+ 116647Wr 6 + 2400r' 

- 2r (m - 2r^j (2IM 7 + 252M 6 r 
+ 1220A4V + 3002M 4 r 3 + 3800.M 3 r 4 
+ 20007W 2 r 5 - 160Mr 6 + 800r 7 ) | , (BIO) 



M (sin 



20r 2 (M + 2r) s 
M 



3M[7M[M+2r) x 



In 



M + 2r 



2M 2 + 3r(cos 29) (SM + lOr) 



+ lLMr + 18r 2 J +6r 2 (cos26»)(^217W b 
+ 259M 5 r + 1288.MV + 3272.MV 
+ 4368.M V + 2704Xr 5 + 320r< 

+ 2r (42M 8 + 609M 7 r + 3778.MV 
+ 12836.M 5 r 3 + 25340X 4 r 4 + 28112.M 3 r 5 
14816A4 2 r 6 + 320CLMr 7 + 1600r 8 ) | . (Bll) 

2. The Magnetic Part of the Weyl Tensor 



c(2) 
're 



3M 2 



20r s {M + 2r) 8 



3840r 8 sin6»cos( 



'(2) 



(sin 26>) [M + 2r 2r 2LM b + 357X 



21707W 4 r 2 + 6342X 3 r 3 + 93887W 2 : 
62167Wr 5 + 720r 6 ) + 21M ( M + 8 



(M + 2r) ]n( 



M 



M+2r 



(B9) 



The magnetic part of the Weyl tensor Bij is defined in 
equation (J2|. To lowest order in the momentum this is 



Ba = B$ + 0(4>) = 4(0)^m i; + 0(4), (B12) 



Jk 



(0)^(1) 



where the extrinsic curvature is given by equation ( B6 ) 



and the covariant derivative and antisymmetric tensor 
are constructed using the unperturbed spatial metric. 
Evaluating equation (B12) results in 



Brd> — Ba 



[M + 2rf 

and the other components vanish to this order. 
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3. Horizon Tendicity 



For the perturbative initial data of Sec. Ill A the ap- 
parent horizon is located at a coordinate distance 



h 



(B14) 



from the center (see equation (24) in [IB])- The inward 
unit normal on the horizon is then given by 



N l = -s 1 = -f 3 Sj, 



where 



X(l,-d e h,0) 



and where the normalization factor A is 

]2 \ 1/2 



A 



IrrleS - {ire)' 



7 rr (<9 e /l) 2 + 2'J r6 dgh + Jgg 



(B15) 
(B16) 

(B17) 



(see, e.g., [H]). We then compute the horizon tendicity 
from 



£ NN = %iVW = %7*V B 



SASB, 



(B18) 



where A and B run over only r and 9. For a diagonal 
spatial metric this simplifies to 



£nn=£ab! aa 1 aa s A a B . (B19) 

This can be expanded in terms of unperturbed and per- 
turbed quantities: 

£nn = R^Xl^+lm-idohfhee) 



+ 7(0)^-7(6)7(0)^ 

+ o(4). 



(B20) 



When evaluated at the horizon location, equation (B20) 
becomes 



£ 



NN 



1 e 2 

k + E nx(- 1711 + 2688 In 2 + 

4.M 2 2560.M 2 v 

3(cos20)(-1871 + 2688 In 2)) + 0(4)-(B21) 

We could plot this function for £nn, but plotting the 
deviation from the proper-area weighted average value of 
£nn over the horizon is more interesting. This average 
is 



cave 
°NN 



IT So E nn^(1 + (d e h) 2 h- 2 ) 1/2 sin 



So So V^ 2 (l + (d e h) 2 h- 2 y/ 2 smOffldt 

(B22) 

which evaluates to 



cave 
C NN 



AM 2 



16M- 



1 



AM; 



> +0(e%), 
(B23) 



as discussed in Sec. IIVC1 

The appearance of the irreducible mass in equation 



( B23 ) is not surprising given that the denominator of 



equation (B22| is 16ttM? t + 0(4)) which follows from 
[31] if we make the usual approximation of replacing the 
proper area of the event horizon with that of the appar- 
cnt h orizon. To see that the the numerator of equation 
(B22) is as simple as —An + 0(e P ), we use the Newman- 
Penrose formalism [32J EH] as discussed in [TJ] to rewrite 
the horizon tendicity as 



■23t(np-\or), 



(B24) 



where ^IZ is the two-dimensional Ricci scalar for the 
horizon, !R means "the real part of" and p, p, A, and a 
are spin coefficients evaluated on the horizon. Equation 
(B24) is defined with respect to a particular null tetrad 
specified in |12j . but the required spin coefficients can 
still be constructed from purely spatial quantities. We 
find that 



1 



/' : 



-TO 



72' 



r D N 



(B25) 



(B26) 



A = — ^^^(Kij-DjNi), (B27) 
V2 



and 



a = — = to* m? (Kij+DjNi), (B28) 
v 2 

where m? — (e J 2 + ie J 3 )/V2 is a leg of the null tetrad con- 
structed from two orthonormal spatial vectors tangent to 
the horizon, and m*° is its complex conjugate. With the 
above expressions, the combination pp — Act turns out to 
be purely real. Moreover, evaluating the spin coefficients 
for a horizon location of the form h = Ai/2 + ef(9) we 
obtain 



pp 



Act 



2048.M 4 



9M 



: (cos 2 #) 



256 



(/-2(cot6W) (f-2d e d e f) 



•0(4) 



(B29) 



Inserting f(8) = — (M/1G) cos 8 from (B14) all the lead 



ing order terms cancel, meaning that, to our order of 
analysis, the integral in the numerator of (B22) reduces 
to the Ricci scalar term 



2tt 



^n^h 2 (i + {d e h) 2 h- 2 ) l i 2 su 



(B30) 

Since the horizon has the topology of a sphere, the Gauss- 
Bonnet theorem guaran tees t hat the integral is 87r, so the 
numerator of equation (B22| is indeed —An + 0(e P ). 
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Appendix C: Binary Black Holes to the presence of hole 2 in addition to the term due to 

the boost of hole 1 

1. Additivity of the Magnetic Part of the Weyl 
Tensor 



We begin by showing that to calculate Bij for the bi- 
nary, we can, to our order of analysis, add the corre- 
sponding single-hole corrections, i.e. 



(CI) 



e 2 P u + 0{ep) = £t+4, ttl (r)+0( e *) 
2s 

= e%(l + Ul (r))+0(e%), (C7) 



The conformal factor is given by equation ( 43 1 , so the 
extrinsic curvature is 



1 2r Cl 2r, 



c 2 



A 



c 2 p 2 

ij 



(C2) 



At a distance of s/2 or more away from hole 2 we can 
use equation (44) to expand and find 



It was shown in [18] that the effect of the boost of hole 2 is 
of higher order. To lowest order, the extrinsic curvature 
near hole 1 is also unaffected by hole 2, as we discussed 
above. We can use the same approach as in Appendix |B| 
to find the nonzero components of the perturbation of £y 
to second order in e p , except that we now replace u with 
1 + u. This calculation yields the nonzero components 



— Ml Ml \ JC1F1 



■0(4,). (C3) 



At these locations K^ 2 2 = + 0(ep), and since we can 
do the analogous analysis at a distance of s/2 or more 
away from hole 1, we can write the extrinsic curvature 
everywhere as 



(C4) 



The background (two unboosted black holes) Christoffel 
symbol is 



pi(0) _ pi(0)Ci _,_ r i(0)C : 



1 F Hujo 2 , n , 2 \ 



(C5) 



where F 



i(0)C 2 
J* 



is 0(ef>) at a distance s/2 or more away 

^(0)Ci 



from hole 2, and rA 1 has the analogous behavior. To- 
gether, equations (C4| and (C5) imply that the covariant 



derivative of the extrinsic curvature is 



E {2) = 



M 



20r 4 (M + 2r) 



3M[ -2r( cos26M(2FM 7 



+ 399X 6 r + 3052X 



5^2 



12194ATr 



4^3 



+ 27344.M V + 33712.MV + 19776A^r 6 
+ 3360r 7 ) - 7M (3 (cos 29) + 1 Vm + 4r 
5 , ^ 



x (M + 6r) (.M + 2r) In ( 



X + 2r 



2r[ 2FM 8 + 399A4 7 r + 3142X e 



+ 12914MV + 29144.M V + 33712AtV 5 
+ 12576.M 2 r 6 - 8160A4r 7 - 1920r 8 ) ) , (C8) 



D k K VJ = DfK^ + D° 2 K?? P2 + 0(ep), (C6) 

where the covariant derivative in the first term on the 
right hand side takes into account metric terms from 
only the first black hole, and the second term is simi- 
lar. From equation ^ the magnetic part of the Weyl 
tensor is therefore equation (CI I, as desired. 



2. The Electric Part of the Weyl Tensor 

Our calculations for near hole 1 are very similar 
to those for the single boosted black hole in Appendix 
[b| Using equation ( 44 1 the conformal factor ( 43 1 can be 



expanded as in equation ( 34 1 , with the difference being 



that the perturbation term now has a contribution due 



c(2) 
-ve 



3M 2 



20r s (M + 2rf 



3840r 8 sin 9 cos 9 



- sin 2(9 J [M+2rj \ 2r [21M G + 357M 5 r 

+ 217CIW V + 63427W 3 r 3 + 9388X 2 r 4 
+ 6216A4r 5 + 720r 6 ) + 21M (m + 8r 



x (M + 2rj 5 In ( 



M 



M + 2rv 



(C9) 
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7(2) 



— - - N8 [ 3M \7M (M + 2rY 
20r 2 (M + 2rf I \ V J 

X ln (j^)( 3 ( COs29 )( M2 + 7Mr 
+ 14r 2 ) + (2r -M)(M + 3r)) 

+ 2r(cos26>) (2I.M 7 + 336M 6 r + 2275M 5 r 2 
+ 8330.M V + 17528M 3 r 4 + 20608A4 2 r 5 
+ 11664Mr 6 + 2400r 7 ^ 

-2r(M- 2r^j (2IM 7 + 252M 6 r 

+ 1180M V + 2602A4V + 2200MV 

- 1200MV - 3360A4r 6 - 480r 7 ) ) , (CIO) 



and 

+ UMr + 18r 2 ) + 6r 2 ( cos 2^ (2IM 6 
+ 259M 5 r + 1288MV + 3272M 3 r 3 
+ 4368A^ 2 r 4 + 2704A4r 5 + 320r 6 ) 

+ 2r (42M 8 + 609M 7 r + 3818MV 
+ 13156.M 5 r 3 + 26140.M 4 r 4 + 28112.M 3 r 5 
+ 11616A4 2 r 6 - 1920Mr 7 - 960r 8 ) ) . (Cll) 
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